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We show that in the hmit of large and positive atom-atom 
scattering length the properties of an atomic-molecular Bose- 
Einstein Condensate (amBEC) are determined by an univer- 
sal energy density functional (EDF). We find that the opti- 
mal conditions for the formation of a stable amBEC are in 
the regime where there are no shallow trimers and the atom- 
dimer scattering length is negative and comparable in mag- 
nitude with the atom-atom scattering length. At tempera- 
tures lower than Tc the chemical potentials for the atoms and 
molecules can be specified independently. Besides three-body 
recombinations processes into dimers of large size, inelastic 
processes involving the formation of deeply bound small size 
molecular states are possible. These inelastic processes do not 
lead to an efficient heating of the amBEC and can be used for 
its mostly non-destructive monitoring. 

PACS numbers: 03.70.+k, 05.30.-d 

The theoretical prediction [1], along with the experi- 
mental confirmation [2] of the fact that one can manip- 
ulate relatively easily the interaction of slow atoms by 
means of external fields, opened the doors to a long series 
of possible experiments. Perhaps among the most spec- 
tacular experimental results obtained lately are the cre- 
ation of the so-called Bosenova [3] and the evidence that 
something related to a mixed atomic-molecular Bose- 
Einstein Condensate (BEC) was created [4], an object 
envisioned by theory several years before [5-7], see also 
Refs. [8,9] for a couple of theoretical analyses of this ex- 
periment. In the present work we identify a new uni- 
versal regime, which appears in the case when the two- 
body scattering length a is positive (a > 0) and large, 
a ^ To, where ro is the effective range, appearing the the 
low energy parameterization of the s-wave phase shift 
kcot6{k) = — 1/a -I- + . . . Typically rg is of the 

order of the so called van der Waals length {Cem/h^y^^ 
[10]. The universality of this new regime is encoded in 
only two parameters, which define the effective Hamil- 
tonian (except for the mass), one of them the two-body 
scattering length a and the other a three-body charac- 
teristic, which can be chosen in various ways. As we will 
show below, in this regime, when na^ <C 1, the proper- 
ties of a cold mixture of atoms and shallow dimers are 
described by a very simple theoretical model. In the other 
limit, when a < and |a| ^ tq, there is another universal 
regime and a new class of dilute quantum liquids exists 
[11]. The basic ingredients of the model are derived from 
the well known Efimov effect [12,13], which pertains to 



the universal properties of the three-body systems with 
large two-body scattering lengths. 

In order to derive the energy density of an atomic- 
molecular gas let us first consider the dilute regime where 
nrQ ^ 1, with n a typical atomic or molecular density 
and a = ©(ro). The details of the potential are not im- 
portant and the atom-atom interactions can be described 
by a contact term 
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where ipl^, ipa are creation and annihilation operators for 
atoms. The value of A2 is determined (after fixing a short 
distance regulator) by the scattering length a. A3 de- 
pends also on a genuinely three-body length scale, de- 
noted below by 03, that can be determined only through 
the value of a three-body observable. Terms with more 
derivatives or higher body forces can be included, but 
their effect on a dilute gas is small as long as the rel- 
evant momenta satisfy the condition fcro ^ 1 [14-17]. 
The Hamiltonian in Eq. (1) contains, in principle, all the 
information necessary to describe the system in the di- 
lute regime, including the molecular states/dimers with 
binding energy 7? jmc? and radius « a when a > 0. In 
the case a ^ rg we are interested in, the presence of 
bound states within the regime of validity of Eq. (1) ev- 
idences the fact that perturbation theory breaks down. 
Since particle-particle diagrams should be summed up 
to all orders, a simple mean-field approximation is not 
legitimate as well and two-particle correlations have to 
accounted for explicitly. Thus, even though Ha contains 
all the information needed to describe molecular states, 
it is inconvenient to use it directly [16]. 

Let us consider now a regime even more dilute, where 
no?' <C 1. The relevant typical momenta are of the order 
p = hk \/na^h/a <C h/a <C h/rg [14]. We will argue 
later that transitions altering the number of molecular 
states are "very slow" in this regime and, consequently, 
for time scales shorter than the transition rate, the num- 
ber of atoms and molecules are separately conserved. An 
appropriate effective Hamiltonian describing the system 
in this regime is 

^2y2 ^2y2 

Ha^ = -i'l^^^a-i^l{^~e2)i^^ (2) 

+ ^Aaa?Al ^V'q + Km^llpLlPatpm + ^ A^^'i/'m^ Vm" 

Both Ha and Ham are applicable in the regime ka <^ 
kro <C 1. Ham has a great advantage over Ha how- 
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ever. Perturbation or mean-field theory is not valid 
for Ha when p ~ h/a, but it is valid for Ham since 
all non-perturbative physics occurring at the scale ~ a 
and leading to the formation of the bound state (dimers 
of size ~ a) and to the Efimov effect (trimers of size 
~ a) is already neatly encapsulated in the constants 



^aay^am and A„ 



To make that more explicit, let 



us determine \aa , ^am. and Xmm in terms of A2 and A3 
by considering zero momentum scattering of two atoms, 
atom-molecule and two molecules. Using Ha these pro- 
cesses are described by an infinite number of diagrams 
shown in Fig. 1, corresponding to the solution of the 
Lippmann-Schwinger, Fadeev and Yakubovski equations 
respectively [12,13,16,18]. The explicit solutions of the 
Lippmann-Schwinger and Fadeev equations are known. 
If we use dimensional regularization [19], where diver- 
gences proportional to powers of the momentum cutoff 
vanish, the same processes are described using Ham sim- 
ply by the Born amplitude, see the Fig. 1. This leads to 
the identification 
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where the values of the universal constants so ~ 1.00624, 
ci « 1.46 and C2 ~ 2.15 are known [20] and ]e2| is the 
true dimer binding energy, ag > is the value of the 
two-body scattering length for which a trimer energy is 
exactly at the atom-dimer threshold and this constant 
is system dependent, see Fig. 2. The specific value 
of 03 > is determined by short-range three-particle 
physics, when all three atoms are in a volume of or- 
der Tq. The functional form of the atom-dimer scatter- 
ing length was established by Efimov [13]. It is natural 
to expect also that the dimer-dimer scattering length 
is of order a, as there are no known four-body shal- 
low states near the 4^2-1-2 threshold and a w Og. 
The value of C3 has not yet been calculated for bosons, 
but for fermions Randeria has shown at mean-field level 
that C3 = 1 [21]. This result, in connection with Grib- 
akin and Flambaum's reasoning [22], leads us to expect 
that in most cases the dimer-dimer scattering length is 
positive. Notice that the matching described above is 
performed in vacuum (length scales of 0{a)) and it is 
not modified by the many-body physics (length scales 
of 0{a/\/ na^) ^ 0{a)). The fact that the dimension- 
ally regulated Ham is perturbative in the fca <C 1 regime, 
unlike Ha, is firmly established [16,18]. 

The leading order contributions to the energy density 
functional are given by the mean field approximation of 
Ham, see Fig. 1 ( NB, a mean-field approximation for Ha 
would be highly inaccurate even at these low densities!). 
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FIG. 1. Matching between the 2-, 3- and 4-particle ampli- 
tudes computed with Ha (left side of the equality) and Ham 
(right side). Only diagrams containing A2-vertices are shown. 
In a description which includes explictly Feshbach molecules, 
each of these vertices stand for a sum of two diagrams, a 
two-atoms^two-atoms diagram and a two-atoms^Feshbach 
molecule^two-atoms diagram. The effective vertices thus 
defined (right side) can then be used to compute the ground 
state interaction energy in the leading order terms in an na^ 
expansion, which is given by the diagrams after the arrows. 



e2f^r; 
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where na,nm are the atom and dimer densities (see be- 
low the discussion on trimers contribution). In terms 
of many-body diagrams built using Eq. (1) this expres- 
sion corresponds to the exact summation of the lin- 
ear and quadratic terms in densities of all the two- 
particle, three-particle and four-particle diagrams, com- 
puted with the initial Hamiltonian Ha- Using well- 
known arguments [10,15], it is straightforward to show 
that higher order corrections, due to mean-field fluctua- 



tions, are controlled by the diluteness parameter • 
1. This EDF expression is remarkable, since, in the Efi- 
mov regime when a ^ rg, it is universal. Only the ex- 
act value of the parameter C3 for bosons is not explicitly 
known. Notice, however, that the condition na^ <C 1 pre- 
cludes its use arbitrarily close to a Feshbach resonance. 
By changing the units of energy to /ma"^ and length 
to a this EDF acquires the following dimensionless and 
thus universal form: 
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with ^ — ci + C2 cot[so ln(a/a3)] tunable. 

In the a > regime only the atom-dimer scattering 
length can change dramatically from negative to positive 
values and even become infinite. If Xaa > 0, Xmm > and 
XaaXmm — X^^ > 0, then a mixed dilute atom-molecule 
system is a metastable gas. This can be seen as well by 
computing the pressure 
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FIG. 2. The energy of an arbitrary Efimov state (thick solid 
line labeled E3) as a function of 1/a. Following Efimov's 
argument we use |i5|^^^sign(E) on the vertical axis instead of 
energy. The two straight lines labeled 3 — > 2 + 1 and 4 — > 2+2 
are the corresponding three- and four-body thresholds. The 
dashed line shows the qualitative behavior of the scattering 
length ttam- 

Phase separation occurs when Aam > V^aa^mm and 
domain wall structures arc expected to appear [23]. A 
(meta)stable trimer phase might also form if the (yet un- 
known) trimer-trimer interaction is repulsive. Another 
type of domain walls can appear in the opposite limit 
[24] . It is easy to see that there is no region where P = 
and thus no stable liquid phase exists, unlike Refs. [6]. 
An amBEC system would be in chemical equilibrium if 

)nm = £2- (9) 

This type of chemical equilibrium can be established 

through the reaction A + A + A ^ A + A2, which has a 
low rate for dilute systems oc ha'^n'^/m.. sec Refs. [25] and 
below. At temperatures T < Tc <x [na-^Y/'^h^ /mo? <C 
[£2! = It? /ma^ , the rhs of the above chemical equilibrimn 
condition is parametrically much larger than the Ihs when 
no? <C 1. The reverse reaction ^2 + ^—*^ + ^ + ^ is 
therefore exponentially suppressed, as it requires an ac- 
tivation energy larger then the thermal energy. Conse- 
quently, these type of collisions thus lead to a very small 
heating (and very slow chemical equilibration of the sys- 
tem) and its ultimate depletion [25]. Since these two 
reactions (and A^ + A2 A2 + A + A, etc. as well) 
are very slow in an amBEC, the chemical potentials for 
atoms and dimers could be specified independently 

Ma ~ ^aa'^a "b ^amf^mj Mm ~ ^am'^a ~t~ ^mm'^m "b ^2* 

In a trap one has to add the trap potentials for atoms 
and dimers in the rhs of these two relations [10]. 

A number of other processes can occur in such a sys- 
tem. It is desirable to consider three body inelastic col- 
lisions involving one or more dimers in the initial state. 
Particularly interesting are however the inelastic atom- 



dimer A + A2 A* + A*^ (see Ref. [26] for a rate cal- 
culation) and dimer-dimer A2 + A2 A* + A* + A^ 
collisions, where A* and A^ stand for a fast atom (with 
its kinetic energy significantly exceeding the thermal en- 
ergy 3T/2) and A\ for a fast dimer respectively (with a 
binding energy of the order 1\? /mr\ ^ h^/ma^). If the 
temperature of the system is significantly lower then the 
dimer binding energy, i.e. T < Tc <^ Tt? /ma?, these in- 
elastic collisions are not leading to a significant heating 
of the system, but mostly to its depletion, in a manner 
similar to the depletion in the BEC of metastable helium 
[27]. If an atom and a dimer collide and this leads to 
the formation of a dimer with a binding energy of or- 
der h? /mr\ ^ fi? /ma? a relative large amount of energy 
is released (NB T < <C 7? /ma^ < /mrl). It is 
conceivable that transitions to various final bound states 
could be resolved and thus the spectrum of fast ejec- 
tiles would be discrete. The momenta of the outgoing 
atom and of the deeply bound dimer are large as well, of 
the order Ti/tq ^ h/a, and such fast atoms and dimers 
will interact relatively weakly with the rest of the atoms 
and dimers. While the slow/thermal atoms and shallow 
dimers interact with each other with cross sections of the 
order of a^, a fast atom and a deeply bound and fast 
dimer will interact with the amBEC with a cross section 
of order Tq ^ a^. One can expect a similar behavior in a 
dimer-dimer inelastic reaction. The fact that the prod- 
ucts of these inelastic processes are fast and weakly inter- 
acting objects can be used for a non-destructive monitor- 
ing of the amBEC. The rates of these inelastic processes 
are obviously controlled by the parameter rg/a"^ ^ 1, 
since in order to form a deeply bound dimer two atoms 
in the shallow dimer have to change their separation from 
distances of order a to distances of order tq. One can- 
not help here but to make a parallel with the similar role 
played by neutrinos (which are also fast and weakly in- 
teracting prodTicts of inelastic processes) in studying the 
interior of the Sun. 

A possible final channel of a dimer-dimer collision 
could be a trimer— atom final state if > 0. If 

the trimer energy is close to the atom dimer threshold, 
then the atom-dimer cross section could be significantly 
larger than the atom-atom and dimer-dimer cross sec- 
tions (jaamj ^ o, amm)- If this is the case, one can expect 
that the atom-dimer elastic collisions would be chiefly re- 
sponsible for the thermalization of the system. However, 
besides the elastic channel there is also an inelastic chan- 
nel present if aam > 0: a dimer-dimer c;ollision leading 
to a shallow trimer and an atom in the final state. This 
process will lead to a relatively small release of energy 
since the trimer energy is very close to the threshold, see 
Fig. 2. With good accuracy the ground state energy of 
a trimer is given in this limit by £3 « £2 — 3?i^/4TOa^„. 
[12] Moreover, the trimer under these conditions is ba- 
sically a very large two-body object of size aam ^ a, 
made of a dimer and a relatively loosely bound atom to 
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it. In this limit the atoms and dimcrs wiU phase sepa- 
rate, since Xam ^ V^aa^mm- A piirc trimer phase could 
appear instead, iff the trimer-triincr interaction is repul- 
sive. When Gam > and Xam < V ^aa^mm dimcr dimer 
collisions with a trimer in the final state will lead to an 
energy release of the order of h^/ma^, and thus to the 
heating of the system, since the momenta of the emerging 
products are of the order h/a in this case. Consequently, 
the best conditions to create a stable amBEC are ex- 
pected when Gam < (no shallow bound trimers and 
thus no A2 + A2 ^ A3 + A reactions) and \aam\ < cl. 
As a last note, there is essentially no overlap between 
present results and those of Refs. [28]. 

Discussions with David B. Kaplan and Lan Yin are 
greatly appreciated. 
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